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The angular momentum mixing in a non-spherical CSC with nonzero azimuthal quantum number 
and therefore violating time reversal invariance has been examined. The mixing is bound to occur 
because of the equal strength of the pairing potential mediated by one-gluon exchange for all partial 
waves to the leading order QCD running coupling constant and the nonlinearity of the gap equation. 
The free energy with mixing is lower than that with p-wave pairing only, but still higher than that 
QQ ' of the spherical CSL state. 
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Q ! I. INTRODUCTION 

It has long been expected that quark matter at high baryon density and low temperature becomes a color 
' (— I '' superconductor(CSC)[l', CSC is characterized by a diquark condensate, which is analogous to the Cooper pair in 
O ^ an ordinary superconductor but the structure of the condensate is much richer because of the nonabelian color and 
I ' flavor charges. 

For very large baryon density, where the masses of u,d and s quarks can be ignored, the ground state is in the 
(— I ' color-flavor-locked(CFL) phase 3]. The situation becomes more involved in moderate density because of the strange 
I— I, quark mass, (3 equilibrium and the charge neutrality conditions, which will induce a substantial Fermi momentum 
mismatch among different quark flavors and thereby reduce the available phase space for Cooper pairing. A number of 
'^l^ ' exotic CSC phases have been proposed in the presence of mismatch, but a consensus point of view of the true ground 
state has not been reached [a 0|- An interesting alternative in this circumstance is the single flavor pairing, which 
^ is obviously free from the Fermi momentum mismatch. Since the attractive interaction of quarks is provided in the 
fS| ■ antisymmetric color-antitriplet channel, the total angular momentum channel must be symmetric in order to insure 
' the overall antisymmetric of the Cooper pair wave function as required by Pauli principle. Therefore, Cooper pair in 
, single flavor pairing should be implemented at a higher total angular momentum, the obvious choice is the p-wave 
■ pairing(J=l), which had been extensively explored in the literatures [sl. lol. [lOl. [TH . 
00 \ The energy gaps in the single flavor pairing can be divided into two categories, spherical gaps (e.g. CSL) and 
• non-spherical ones (e.g. polar phase and A phase). While a spherical gap is made of p-wave alone, a non-spherical 
gap may contain higher partial waves because of the nonlinearity of the gap equation. It was argued in the literature 
that the contribution of higher partial waves is of higher order in g with g the running coupling constant of QCD. 
This, however, is not the case, as we shall explain. 

In a previous paperfT^, we considered the single flavor CSC with longitudinal pairing, in which the pairing two 
quarks have equal helicity. Because of the equal strength of the pairing potential mediated by one-gluon exchange for 
all partial waves to the leading order QCD running coupling constant, a non-spherical pairing receives contributions 
from all odd J's at the same order of g. A consistent solution to the gap equation with a deflnite azimuthal quantum 
number A/, however, can be constructed. For the CSC that pairs red and green quarks only, the energy gap with a 
definite M is of the form 

A = e'*^'^Ao/(cos0) (1) 
with 9 and ip the polar angles of the relative momentum p of the pairing quarks, where 

/(cos0)= bjP^\cos9) (2) 

J>|AJ"|,J=odd 
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with Pj^(x) the associate Legendre function. The free energy will be brought down by the mixing in comparison with 
that reported in the literature. For M — 0, we found that 

/(cos 61) = e""/^[0.9866Pi(cos6i) - 0.0673P3(cos6i) + 0.0214P5(cos6i) + ...]. (3) 

The magnitude of the condensation energy is raised by 3.5 percent, which is too small to compete with the spherical 
color-spin-locked(CSL) state in single CSC. As we pointed out in [ij], the energy balance between the CSL and non- 
spherical states maybe offset by the presence of anisotropy or the violation of time reversal invariance. In a compact 
star, the magnetic field and the stellar revolution provides such an opportunity. 

In this paper, we extend the analysis in to case with a nonzero azimuthal quantum number (A/ ^ 0). We 
will firstly derive the nonlinear integral equation for the angular dependence of the gap function and then give the 
numerical solution to this equation. Finally, a summary and concluding remarks will be given. 



II. THE INTEGRAL EQUATION OF THE ANGULAR DEPENDENCE 

The non-zero azimuthal quantum number violate the time reversal invariance and the gap function is proportional 
to e^"^ in this case[14j. The condensation energy density for single flavor CSC with longitudinal pairing reads [T^ 



F= - ^ jdy j dy' j d'p J d^p'V{y~y\Q) x 



dv d'p - ,,,,,, (4) 



(27r)3y J ^v^ + \4>{i,,p)\ 

where v and v' are Euclidean energies, Q is the angle between p and p' and V is the pairing potential mediated by 
the one-gluon-exchange of QCD 

V{v - v\ 6) = Di{v- v\ 6) + Dt{v - u' , 9) (5) 

with Di and Dt is the longitudinal(color electric) and transverse (color magnetic) part of hard-dense-loop(IIDL) gluon 
propagator respectively. The gap function p) is extracted from the Nambu-Gorkov off diagonal blocks of the quark 
self-energy and its value at u = Q gives A of ([1]). The gap equation can be derived by minimizing the condensation 
energy with respect to the gap function 

and we end up with 

A consistent derivation of the gap equation up to the subleading order need to include the one-loop self energy of 
quarks, the net result is to replace the first term in the square root on RHS of Eq.Q by v'"^ jZ'^iy'^ with Z{y'^ the 
wave function renormalization factor (islllGj. But it will not interfere with the angular dependence of the gap function 
as we have seen in[l2|. 

The potential V can be expanded in terms of Legendre polynomials [lB| 

V{v~v',Q) = — ln--^^(2J+l)Pj(cose) + ^^(2J+l)cjPj(cose) (8) 
with ujc = ^"^^^^''^ and cj given by 



J 



2E-- (9) 

z — J ri 



1 

n— 1 



From this expansion, we can see that the pairing strength are equal to leading order of the QCD running coupling 
constant, but subleading terms fall off with increasing J. It is this falling off that makes the amount of the angular 
momentum mixing numerically small for the solutions in(l^ and in this letter as we shall see below. 



Substituting Eq. ([8]) into the gap equation ([T]), we have 



21: J 1 ~ p ■ p' 

where — /{I^tt"^) and a UV cutoff loq ^ gfi is introduced. On writing 

^(z/,p) = e**^'^V(i^,cos6'), 
the equation for the real function ipi'^: ^) reads 



f /d.'{i(ln^^+ln^) 



3 y dx' 



where x = cos and 



2^ /.27r ^iM(ip'-^p) 

Im{x,x') ^ tt I ' — ~ = '■^'^1 



27r 



1 — p • p' 



The most favored pairing channel corresponds to \M\ = 1 and there is no difference between M and 
consider the case of M — 1 below. We have Ii (x, x') = I{x, x')/\x ~ x'\ with X given by 



1 — xx' — \x — x'\ 



^(1 -a;2)(i _a;'2) 

The gap equation (|12p can be further simplified by using the approximation in[ 



hi I r — hi- 



with \v-^\ — max(|:/|, Introducing new variables 



4 = In — , a = in — 



and writing ip{v,x) as we have 



|a; — x'l 



where $(^, x) is defined by 



/•OO 



^iC,x) 



^1 + q|^e2C' 

Taking the first derivative of both sides with respect to ^ in Eq. (fT7|) , we have 



which implies 



4 



as ^ ^ oo for all x. Another derivative of ((T9|) yields 



0. 



The x-dependence of '4'{£,,x) will be fixed by the Eq.([T7]) at — a, i. e. 

„ /"^ , ,^(a,x')I - ^(a,x) 
a^(a,x)-^(a,x)+3 dx ^ ' / ' = 



(21) 



(22) 



The solution to ((2T|) subject to the boundary condition ((20)) has been obtain inllj] up to subleading order, which 
applies to the case Af = 1 as well. We have 



Tl;{a, x) ~ Ao/(x) [- - 5(6 - a) ~ g\n2\ + 0(5) 



and 



where Aq is the s-wave gap given by 



$(a,x) =5Ao/(x) + 0(g) 



2 ^ Ao 



(23) 



(24) 



(25) 



and h = In 
dependent factor f{x) 



Substituting Eg. ([^5]) and Eq. ([M)) into Eg. ([2^ . we derive the integral equation for the angle 



/(.)ln|/(x)|+3/\x- ^(-)"/^-')f '-') ^0 

-1 \X ~ X \ 



(26) 



III. THE SOLUTION TO THE INTEGRAL EQUATION 

As we have done infT5], this type of integral equation can be solved by a variational method. Substituting the 
solution of ([2T|) and ([20)1 into (|4]), the condensation energy density becomes a functional of /(x)(see appendix B 
indl), 



F = 



27r2 



with 



m{x)]^ I ^dxf'{x)[\n\f{x)\-- 



dx / dx 



,P{ x)-2f{x)f{x')I+p{x') 
\x — x'\ 



(27) 



(28) 



It can be easily verified that the variational minimization of Eg. (P5|) does solve Eq. ipS)) . Before the numerical calcu- 
lations, we consider a trial function as 



fix) = cP^{x) = cy/l - x^ (29) 
Substituting it into the condensation energy density Eg. (|28p and the minimization with respect to c yields 

1 



(30) 



at which 



-5.863 X lo- 



cal) 



This trial function is what people carried over from A phase in '^He[T3, [HI, which contains p-wave only, but it is not 
optimal. The free energy will be lowered further by including higher partial waves as we shall see. 
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FIG. 1: The angular dependence of the gap function with nonzero azimuthal quantum number. The dashed line and solid one 
correspond to the trial function and the numerical solution to Ea. (|26p respectively. 

Following the same procedure infl^, we obtain the numerical solutions to Eq. ipS)) . In FiglH we show the angular 
dependence of the gap function with nonzero azimuthal quantum number, the dashed line and the solid line are the 
trial function and the numerical solution to Eq. (j26p respectively. They depart from each other slightly, indicating a 
small mixture of higher partial waves. We find the minimum value of the target functional with the solid line 

= -5.977 X 10"^ (32) 

which drop from Eg. (1311) bv 1.9 percent. 

As we have found in|12|. the drop of the condensation energy with longitudinal pairing with zero azimuthal quantum 
number by angular momentum mixing is numerically small. Here, with nonzero azimuthal quantum number, the 
situation is also true. With these small drop amount in condensation energy, the non-spherical pairing in single flavor 
CSC can not compete with the spherical pairing state CSL. 

Regarding the angular momentum contents for our solutions according to Eq.Q for M = 1, we found 

/(cos6') = ie^/^e"^[1.0112Pii(cos6')+0.0212P3i(cos6i) + 0.004^5^(0086') + ...]. (33) 

IV. CONCLUDING REMARKS 

In summary, we have explored the angular momentum mixing in a non-spherical CSC with nonzero azimuthal 
quantum number, which violates the time reversal invariance because the gap function is not real in this case. The 
mixing is driven by the equal strength of the pairing potential mediated by one-gluon exchange for all partial waves 
to the leading order QCD running coupling constant and the nonlincarity of the gap equation. However, the gaining 
factor of the condensation energy caused by mixing is smaller than that by forming the spherical CSL state. The 
angular momentum mixing in the parallel case with transverse pairing has also been considered in 13], in which the 
amount of the free energy drop is also small. Therefore, we conjecture that angular momentum mixing of various 
non-spherical CSC is not sufficient to compete with the CSL energetically in the ultra relativistic limit. A rigorous 
proof of this statement is anticipating. 
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